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Abstract 

We study the finite-size behaviour of a tree-level on-shell improved action for 
the iV-vector model. We present numerical results for N = 3 and analytic results in 
the large-iV limit for the mass gap. We also report a perturbative computation at 
one loop of the mass gap for states of spatial momentum p. We present a detailed 
comparison of the behaviour of this action with that of other formulations, including 
the perfect action, and a critical discussion of the different approaches to the problem 
of action improvement. 



1 Introduction 



Lattice simulations are at present the most effective method to investigate non-perturbative 
properties of field theories like QCD. By necessity Monte Carlo simulations are performed 
on finite lattices and at finite values of the correlation length. It is therefore of utmost 
importance to understand the systematic effects due to scaling corrections. 

In the last years a lot of work has been devoted to invent lattice models that have 
small scaling corrections so that continuum results can be obtained on small lattices and 
thus with a limited use of computer time. 

This program was started by Symanzik |3J who put up a systematic method to improve 
asymptotically free theories using perturbation theory. Soon after, Liischer and Weisz 
noticed that simpler actions could be used if one was only interested in on-shell quantities 
. A different approach is the perfect- action program in which improved actions are 
determined as fixed points of renormalization-group transformations |3|-|9|]. By definition 
classically perfect actions do not show any lattice effect at tree level. In the standard 
language they are Symanzik tree-level on-shell improved to all orders of a Jl0[ . 



Recently the improvement program has been implemented non-perturbatively for the 
fermionic action [rT-ITjj]. This represents an important step forward. Indeed actions that 
are improved to a finite number of loops have scaling corrections of order 0(alog p a), or, 
in the statistical mechanics language, of order 0(£ _1 log p £), where £ is the correlation 
length. On the other hand non-perturbatively improved actions should have corrections 
of order 0(£~ 2 log 9 £). It must be noted however that improvement comes with a price. 
Improved actions are more complicated than unimproved ones. Thus to understand the 
practical relevance of any improvement one should also consider the additional cost in the 
simulations. 

In this paper we will study the two-dimensional iV-vector model. This theory pro- 
vides the simplest example for the realization of a nonabelian global symmetry. Its two- 
dimensional version has been extensively studied because it shares with four- dimensional 
gauge theories the property of being asymptotically free in the weak-coupling perturbative 
expansion |l5|-p!7[| . This picture predicts a nonperturbative generation of a mass gap that 
controls the exponential decay of the correlation functions at large distance. 

Besides perturbation theory, the two-dimensional iV-vector model can be studied using 
different techniques. It can be solved in the N = oo limit [IE, |T^] and 1/N corrections 
can be systematically calculated 



20-22]. An exact S*-matrix can be computed [23,24" 



and 
limit (3 



using the thermodynamic Bethe ansatz 

11,11 



oo has been obtained 



the exact mass gap of the theory in the 
The model has also been the object of extensive 



numerical work p7-32| mainly devoted to checking the correctness of the perturbative 



37 . It is on-shell tree- 



predictions fl33|-|36[ . 

We will consider the action that has been proposed in Ref. 
level improved and satisfies reflection positivity so that a positive transfer matrix can be 
defined. Moreover we will show here that this action can be efficiently simulated: indeed 
one can use a Wolff algorithm |3"E|-pET| with unfrustrated embedded Ising model. Therefore 
no critical slowing down is expected. This is at variance with the standard Symanzik tree- 
level action. Indeed, also in this case one can define a cluster algorithm [42]. However the 



embedded Ising model is frustrated and therefore critical slowing down is still present. 
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The purpose of this work is twofold. First of all we want to understand quantitatively 
the effect of tree- level improvement. Indeed it is not obvious a priori that the idea is 
effective since the corrections to scaling are simply reduced a logarithm of a. As we 
shall see from our numerical results for N = 3 and our analytic expressions in the large- 
iV limit, tree-level improvement effectively works: indeed one finds "naturally" - - i.e. 
without any additional tuning of the parameters — a reduction of the scaling corrections 
by approximately a factor of two. The second purpose of this work was to understand the 
results of Ref. for the perfect action which shows, for a particular value of the parameter 
k that parametrizes the renormalization-group transformations, a dramatic improvement 
with respect to the standard action. The on-shell action that is considered in this paper 
differs from the action of Ref. || by terms of order a 4 and it shares the property of being 
extremely local. We thus wanted to understand if the results of Ref. |J depended only on 
the improvement and on the locality of the action. The answer is clearly negative, since 
the on-shell action we study has much larger corrections. At this point the question that 
arises naturally is whether the exceptionally good behaviour is related to the fact that 
the action is a classical fixed point of a renormalization-group transformation. Extending 



the perturbative calculation of Ref. |10[ we will see that this interpretation is unlikely. 
Indeed with a different choice of the parameter k one can define perfect actions that are 
local but that are not expected to have such a good behaviour. 

In this paper we will investigate the corrections to finite-size scaling (FSS) for various 
models that are introduced in Sect. |2|. A discussion of the FSS corrections in the large- N 
limit is presented in Sect. |[ We will show here that tree- level improvement reduces the 
corrections to scaling from log L/L 2 to 1/L 2 . Sect. [| presents our algorithm and Monte 
Carlo results. Finally in Sect. [|we present our conclusions and compare our results with 
those obtained with other types of action. App. |A] presents the details of the analytic 
calculation of the FSS functions for the on-shell action in the large- N limit, App. [B| the 
analytic computation at one loop of the mass gap for non-vanishing spatial momentum 
while in App. |C] we give some details of our calculations with the perfect action. 



Preliminary results of this work have been presented at the Lattice 97 conference [43 



2 The models 



In this work we will study in detail the FSS properties of the tree-level on-shell improved 
action proposed in Ref. pTj : 
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(2.1) 



where d runs over the diagonal vectors (1, ±1), while in the last term 1 and 2 are the unit 
vectors along the x- and y-axis respectively. It was shown in Ref. |J7| that the action 
( |2.1| ) is reflection-positive and tree-level on-shell improved. In the formal continuum limit 
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we have 



S 



onshell , 



d 2 x 



2 V ; 2 



(<r ■ DV - (<r ■ Dcr) 2 ) + 0(a 4 ) 



(2.2) 



where □ = ^„c^ and a 2 = — ^. 

We will study the FSS corrections for this action and we will compare these results 
with those obtained for other actions that have been extensively studied in the literature: 



the standard action 



S 



stdi 



<j x ■ cr x+fl 



(2.3) 



whose FSS behaviour has been extensively studied in Refs. ]2U|. [Ti|. |"Ti 

the Symanzik action M 
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X,fl 
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-<J X ■ <J x +2fj, 



(2.4) 



which is designed to cancel 0(a 2 ) lattice artifacts in tree-level Green functions [| 



the diagonal action [46 



S dla9 (cr) = ^2 f ■ Ca;+/i + ■ CTx+d j 



(2.5) 



whose one-particle spectrum has no 0(a 2 ) artifacts at tree level but that is not on- 
shell improved. For instance the four-point function shows 0(a 2 ) corrections even 
on- shell. 



the classically perfect action [g|,|J that is defined as the fixed point of a class of 
renormalizat ion-group transformations for (3 = oo. The two-spin part is given by 



(2.6) 



where w(x) is the so-called perfect laplacian ||7|- Its Fourier transform is given by 



w(q) 



(gi + 27r/!) 2 + (g 2 + 2vrZ 2 ) 2 ( ?1 + 27r/i) 2 (g 2 + 2vr/ 2 ) 2 



(2.7) 



where k is a parameter that characterizes the renormalization-group transformation. 
The four-spin coupling — and also higher order couplings — cannot be computed in 
closed form. In App. |C| we give some details on our determination of the four-spin 
term for various values of k. The formal continuum limit of the perfect action is 
given by Eq. ( |2.2|) with a 2 = (k, — 4)/(12k). Notice that the perfect action with 
k = 2 and the on-shell action (|2.1| ) are, at tree level, equivalent also off-shell up to 
terms of order a 4 . 
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3 FSS functions in the large- N limit 



In this section we will discuss the relation between improvement and finite-size scaling in 
the large- N limit extending the discussion of Ref . [ 1HJ . We will consider anlxT square 



lattice and we will generalize the action (|2.1| ) considering 

S[cr] = N^J(x - y)cr x ■ cr y + ^ K(x - y)K(x - z)(cr x ■ a- y )(cr x ■ cr z ), (3.1) 

x,y x,y,z 

where 

K(p) = J2 K ( x ) e ~ ipX =P 2 + °(P 4 )' ( 3 - 2 ) 

X 

J(p) = J2j^)e- ipx = J(0)-\w{p), (3.3) 

X 

w(p) = p 2 + «!^^ + a 2 (p 2 ) 2 + 0(p 6 ). (3.4) 

Here p 2 = J2^Pfi = sin(p A1 /2)) 2 . 

We will assume, as usual, that all the couplings respect the symmetry of the lattice and 
that w(p) vanishes only for p = in the Brillouin zone. Moreover we require the action 
to be ferromagnetic, that is to have a unique maximum corresponding to the ordered 
configuration. The general class of actions (|3.1| ) was studied in Ref. |37j where it was 



shown that S[cr] is tree-level on-shell improved if a>i = 1/12 and a2 = a^. 

In the large-iV limit the model can be solved using a standard Lagrange-multiplier 
technique. One introduces two parameters mi,T and u>l,t related to (3 by the gap equa- 
tions)]: 

P(l + u L , T ) = TZ — 2 — ' (3-5) 

P hM \ 1 sr *M fie* 

- — {1 + u L: t)uj L! t = T7^Z^~( ZZ — 2~' 3 - 6 

2a 3 LT ^— ' w(p; u LiT ) + m l L T 



where 



_ w{p) + ujK(p) 

w{P, = — • 3.7) 

l+oo 



The two-point function is simply given by: 

1 1 gipC^-a) 



(J,,) = (3(l + u L , T )LT\w(p;u L , T )+ml T (i 



•y 

>y i ■ ■ /' i u i ■ 

P 



1 We should remark that for some choices of K(x) and J(x) these equations may not be correct for 
all values of j3 (see for instance the exact solution of the mixed O(N) — RP N ^ Y models of Ref. p9[). 
However under our hypotheses they should always describe the theory for sufficiently large values of (3. 
This is sufficient for our analysis since we are only interested in the critical limit /3 — > oo. 
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In Appendix [A] we present an analytic computation of the FSS curves and of their leading 
correction. For the mass gap [i(L, (3) defined on a strip L x oo, in the limit (3 — > oo, 
L — > oo, with /x(L, (3)L = x fixed, we find 



a„(s;L) 

L 2 



+ 0(L~ 4 \og 2 L) 



(3.9) 



The function is the FSS function and it was already computed by Liischer f50| 



f^x) = exp 
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71=1 
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exp 
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;i + o(x)). 



(3.10) 



The function A M (x;L) is the leading correction to FSS. For L — > oo keeping a; fixed, it 
behaves as log L and it has a regular expansion in powers of 1/ log L of the form 



A M (x; L) 



oo 



x)(logL)' 



(3.11) 



The leading term is given by 



S n,-i( x ) 



x 

T 



;i-/ M (x))(l-12a 1 -16a 2 + 16a 3 



(3.12) 



It vanishes for tree-level improved actions. Thus improvement reduces the corrections to 
scaling from log L/L 2 to 1/L 2 . Generically, the next coefficient S^o (x) does not vanish 
for tree-level improved actions. However one can choose K(x) and J(x) in order to have 
fifj.,o(%) = (see Appendix |A|). In other words, with an appropriate choice of these two 
functions it is possible to improve the action to one- loop order. One can verify that it 
is not possible to obtain <^,i(a;) = for a 3 ^ 0. In this case, actions of the form (ft.ip 
cannot be two-loop on-shell improved. This is not unexpected since, six-spin couplings 
are needed for two-loop improvement. If 0:3 = the situation is simpler. Indeed in this 
case, once the action is one-loop improved, it is automatically improved to all orders of 
perturbation theory |37]. This is however an accident of the large- N limit and it is not 
true for finite values of N. 

Since tree-level improvement reduces FSS corrections only by a logarithm of L, one 
could be skeptical on the effectiveness of the idea. We have therefore compared the FSS 
scaling corrections for various actions. In Table [l] we report /i(oo, (3)/fj,(L, f3) for various 
values of L for x = 2 for the on-shell action ( |2.1| ), as well as for the actions ( |2.3| ) and 
(|2.4j ). It is clear that tree-level improved actions show much smaller corrections. For the 
purpose of future comparison we consider also the quantity 



R(L,x) = 2L/i(2L,/3). 



(3.13) 



It is easy to convince oneself that in the FSS limit R(L, x) assumes a finite value R(oo, x) 
with corrections of order logL/L 2 . In Fig. [I] we show R(L,x) for x = 1.0595 for various 
actions. Also here it is evident that improved actions show smaller corrections to scaling. 
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Cjstd 


Cjsym 


7Z r — n — 

oonsnell 


4 


0.79362345 


0.78565020 


0.78324283 


6 


0.78640181 


0.78100623 


0.78024309 


8 


0.78290367 


0.77913556 


0.77884796 


10 


0.78099463 


0.77822405 


0.77808048 


16 


0.77857473 


0.77720115 


0.77717924 


20 


0.77792482 


0.77696043 


0.77697663 


32 


0.77714433 


0.77669433 


0.77669910 


oo 


0.77652331 


0.77652331 


0.77652331 



Table 1: /x(oo, /3)//i(L, (3) for x = 2 for different lattice actions. 



Let us now discuss the validity of the expansion (|3.9|) . First of all the expansion is not 
uniform in x: the error increases as x — > oo. This fact can be checked analytically from 
the exact expressions of App. |A]. Its origin is easy to understand. In order to have FSS 
one should work at values of (3 such that the correlation length is much larger than one 
lattice spacing, i.e. fi(L) <C 1. This means that the expansion (|3.9| ) makes sense only for 
x/L < 1. 

We want finally to discuss the connection between the FSS limit for x — > and 
perturbation theory (PT). The PT limit corresponds to (3 — » oo with L fixed. For instance 
for fi(L) one obtains 

xsl tL)L=±£?gL, (3.14) 

^ n=0 " 

where a n (L) has an expansion of the form 

n 1 n 

<L) = £ affi log fc L + «2! ^ L + 0(L' 4 ). (3.15) 

k=0 k=0 

This expansion can be inverted to give 

1 oo 

(3 = -J2b n (L)x n , (3.16) 



X 

n=0 



where the coefficient b n (L) have an expansion of the form ( |3.15 ). To obtain the FSS curve 
one can then use the asymptotic freedom prediction 

and substitute the expansion of (3 in terms of x. In this way one computes /i(oo)//x(L). 
The leading term for L — > oo is L-independent and correctly reproduces the expansion of 
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the exact result ( p. 10 ). One can also compute the first correction in powers of 1/L 2 . In 



this way one obtains the expansion of A At (x; L) for L fixed in the limit x — > 0: 

A^x-L) = ^P n (logL)x™, (3.18) 

n 

where P n (x) is a degree-n polynomial. Clearly this expansion is incorrect for L — > oo, 
although, for finite small x, it gives a reasonable approximation as long as L <C L max (x) ~ 
e 7r/a;_ j n ^ e pp 0S it e limit L ^> L max (x), the correct expansion is given by Eq. ( |3.11| ). 
This result reflects the fact that the PT limit (3 — > oo at fixed L (equivalent to x — > at 
fixed L), followed by L — > oo does not commute with the FSS limit followed by x — > 0, 



except for the leading term [48 



Numerically PT calculations of the corrections to FSS can still be useful, at least if one 
is able to compute enough terms. We have considered for instance R(L, x) that admits 
an expansion of the form 

R(L,x) = Rq(x) + —,R x {x]\ogL) + . . . (3.19) 

where Rq{x) = R(oo,x) is the FSS function. For the on-shell action (|2.1|) we have up to 
three-loops 

R x [x\ logL) = -0.0251916x 3 + (0.0321612 + 0.008018747 logL)x 4 + 0(x 5 ). (3.20) 

Because the action is improved there is no x 2 (one-loop) term. Consider now AR(L, x) = 
R(L, x) — Ro(x). For x = 1.0595 and L = 8, 64 we obtain the estimates 

AR(8, 1.0595) = (1 loop); -0.4681 • 10~ 3 (2 loops); 0.4934 • 10~ 3 (3 loops), 
Ai?(64, 1.0595) = (1 loop); -0.7314 • 10~ 5 (2 loops); 0.1284 • 10~ 4 (3 loops). 

This should be compared with the exact results 0.751 • 10~ 3 for L = 8 and 0.147 ■ 10~ 4 for 
L = 64. The three- loop result is reasonably close to the exact value while the two- loop 
approximation gives a grossly inexact guess (it has even the wrong sign) of the exact 
result. 



4 Monte Carlo results 



As we have seen in the previous section tree-level improvement changes the correction 
to FSS by a logL. A priori this appears to be a very small improvement. However at 
N = oo the actions (|2.1| ) and (|2.4j) show a much better behaviour with respect to the 
standard action ( |2.3|) . We decided to check if this is true also for N = 3, comparing our 
data for the mass gap with the simulation results of Ref. fHJ for the standard action and 
of Ref. for the perfect action. 



We have used in our simulation the Wolff algorithm p8|-|41|| with standard Swendsen- 
Wang updatings [^T],|52[. The idea is the following: given a spin configuration {(t x }, one 
chooses randomly a unit vector r and defines Ising variables e x rewriting 



(4.1) 
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where tr\\ x and <r± x are the components of er x respectively parallel and perpendicular 

to r. Eq. ([4. ID defines an effective action S e ff({e}) for the embedded Ising variables. 

A simplifying feature of our action fl2.1|) is that the effective action ^•"({e}) contains 
only two-spin Ising couplings and therefore it can be updated using a standard cluster 
algorithm. We can indeed rewrite the action ( |2.1| ) as 



S 



onshell 
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The Wolff embedding acts on these couplings as follows: 
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{(T X ■ (Ty) 2 — > 2e x €y(T\\ x ■ CT\\yCT\_ x ■ (T±y 



The embedded action becomes 

S 6ff ({e}) = £ 



where d are the diagonal vectors (1, ±1). The couplings are defined by 



J (1) 

"XfJ, 



1 



4 1 

2 ~~ 3°"^ ' <T±X +^ 



(4.2) 
(4.3) 



{(T y ■ cr x )((Ty ■ a z ) — > e x e z tT\\ x ■ (T\\ y (T\\y • <X|| Z (4.4) 

+ ej/CzCTllj/ ' C\\z (T ±x ■ C ±y 

+ tyt X CT\\y ■ <7\\ X <7 ± z ■ O ±y 



(4.5) 



12 



cr±x ■ cr± 



J 



(2) 



1 



6 



"II* 



12 

cr± x ■ cr ± x -v+n 
1 , 



12 



& ±x ■ & ±x-u~ 



1 — Cm , 1 — , , 1 

2 If+i 2 ll*+« — 1 



(4.6) 
(4.7) 



where, in the first definition, v is a unit vector orthogonal to fi. It is easy to see that the 
embedded Ising model is not frustrated. Indeed if we redefine the Ising variables as 



e' x = e x sign (cr • r) (4. 

rW ^ n 7(2) 



we obtain a new effective action of the type (|4.5|) with Ji M > 0, J^J > 0. Therefore we 



expect the algorithm not to show any critical slowing down [41 
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The purpose of our simulation was the determination of the FSS function R(L, x) 
defined in Eq. ( |3.13|) for x = //(L, (3)L = 1.0595 for a strip L x oo with periodic boundary- 



conditions (pbc) in the spatial direction. This is indeed the value of x for which results are 
available for the standard |44| and for the perfect action ||. Given L, in order to compute 



R(L,x), we have first determined f3 such that fi((3,L) = 1.0595. This has been obtained 
performing a Monte Carlo simulation at a nearby value f3 run of /3, and appropriately 
reweighting the correlation functions. Analogously f/,(/3, 2L) has been determined from 
a Monte Carlo simulation at a nearby value of (3 and then applying the appropriate 
reweighting. Since we use small lattices, the reweighting technique works very well and it 
does not increase significantly the error bars. 

Since we wanted to determine the mass gap with high precision we paid particular 
attention to systematic effects. Most of the simulations were done on lattices L x T with 
pbc and T = 10L. We computed the correlation function 



G(t) = J^J2 ■ v^v+t), (4-9) 

x 1 ,x 2 ,y 

and extracted the mass gap assuming 



G(t) ~ coshm - tj (4.10) 

for L < t < T/2. To check for possible systematic deviations we also performed two 
simulations with "cold wall" boundary conditions: in this case the spins at one temporal 
boundary (t = 0) where fixed in one direction cr waU while we used free boundary conditions 
on the other side. Of course, in the spatial direction we still used pbc. In this case the 
mass gap was extracted from 

G(t) = J2(<r waU ■ *t*) (4-11) 

x 

that is expected to behave as e~ mt for t large. This type of boundary conditions auto- 
matically projects out the excited modes that have energy 

Ehtf, L) - Eotf, L) = ^jlil + 1) + 0(/T 2 ), (4.12) 

and therefore reduces the systematic errors. However it has a disadvantage for our pur- 
poses: with respect to pbc, much longer simulations are required to obtain the mass gap 
with the same precision. We have performed two runs at L — 5, 10 with "cold wall" 
boundary conditions and we have not observed any systematic difference with respect to 
the runs with pbc (see Table |2|). Therefore we believe that our results have a systematic 
error that is smaller than the statistical one. 

The results of our simulations are reported in Table |2|. The total simulation took 
approximately 5 months of CPU-time on a SGI Origin2000. 

These results are shown in Fig. |2] together with the previously obtained ones referring 



to the standard ||4| and perfect actions. Clearly the on-shell action is better than the 
standard action. However its behaviour is worse than that of the perfect action which 
shows no scaling corrections even for L = 5. 
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LxT 


b.c. 


Prun 


a 


lW)L 


N s tat 


Tcpij(ms) 


5 x 30 


periodic 


1.212 


1.21045 ±0.00025 


* 


15 • 10 B 


11/4.5 


10 x 100 


periodic 


1.212 




1.2780 ±0.0019 


3.75 • 10 6 


78/30 


5 x 50 


cold wall 


1.212 


1.2104 ±0.0004 


* 


84 • 10 s 


11/8 


10 x 60 


cold wall 


1.212 


* 


1.278 ±0.002 


175 • 10 6 


25/18 


7 x 70 


periodic 


1.260 


1.2691 ±0.0004 


* 


21.2- 10 e 


38/15 


14 x 140 


periodic 


1.260 


* 


1.2782 ±0.0009 


11 • 10 6 


180/60 


10 x 100 


periodic 


1.340 


1.32835 ±0.00018 


* 


16.2- 10 B 


78/30 


20 x 200 


periodic 


1.330 


* 


1.2738 ±0.0007 


7.2 • 10 6 


420/122 



Table 2: Monte Carlo results. (3 run is the value of (3 at which the Monte Carlo simulation 
was done, $ is the value of f3 such that L) = 1.0595. N stat is the number of iterations 
of a Wolff algorithm with standard Swendsen-Wang updatings. Tcpu is the CPU time in 
ms for a iteration on a SGI Origin2000: the first number is the total CPU time spent in 
our simulation, the second number is the time spent in the updating only, i.e. without 
measuring the two-point function. 



5 Critical discussion 

In this paper we have investigated the FSS behaviour of the action fl2.1|). The main 



point was to understand, in a model in which high-statistics data can be generated, 
how effective perturbative improvement is. Our large- N analysis presented in Sect. ^ 
shows that tree-level improvement effectively reduces the scaling corrections. In this 
limit, simulations with the standard action on a lattice of volume L 2 would be affected by 
the same systematic errors of results obtained on a lattice of size of order L/2 using the 
action Q (cf. Table | and Fig. 0). For N = 3 our numerical data show a similar effect. 



Therefore, using the tree-level improved action ( |2.1| ) one can perform simulations on a 
lattice that is four times smaller and still obtain the same scaling corrections. Of course, 
in order to make a fair comparison of the two actions, we should take into account the 
fact that the action ( |2.1|) is more complicated and therefore one Monte Carlo iteration is 
slower. For instance, one update with the on-shell action takes 30 ms on a 10 x 100 lattice 
(see Table |2|), while it takes 11 ms if one uses the standard action. Therefore, the time 
spent in the updating is reduced only by 20-30%. However it is also important to take 
into account the time spent in the measurement of the observables. In our simulation we 
measured the two-point function: the total CPU-time per iteration on 10 x 100 lattice is 
78 ms, which should be compared to 300 ms spent by the standard action on a 20 x 200 
lattice. Using the on-shell action we obtain comparable results for this quantity in 1/4 
of the CPU time. Notice that we expect the two-dimensional cr-model to be the case 
in which Symanzik improvement has the smaller pay-off. Indeed in this case the lattice 
volume increases only as L 2 , and the algorithm does not have critical slowing down, so that 
working with large volumes is not so expensive. In QCD the situation would be radically 
different since in this case the CPU time to produce an independent configuration scales 
as L A+Z with z >2. Thus, even a small reduction of the needed values of L, significantly 
reduces the computer time. 

Of course, the skeptical reader may think that the better behaviour of the improved 
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action is a mere numerical coincidence. After all, we have simply improved the behaviour 
of the corrections to scaling by logL, which is a very slowly varying function: therefore, it 
is not even obvious that, for our range of values of L, say L m 16-128, an improved action 
behaves better than the standard one. Indeed it is possible to write down (complicated) 
actions that are not tree-level improved and yet, for L « 16-128, behave better than the 
action ( [2.1|) , and, viceversa, to invent improved actions that, in the same range of values 
of L, are worse than the standard action. However the relevant question — unfortunately 
a not well defined one — is whether "simple" local tree-level improved actions behave 
better as long as L is larger than 1. In order to understand if there is a positive answer, 
we should study different actions. We have considered the Symanzik tree-level action 



2j4j ). For N = oo the actions (|2.1|) and (2.4) behave similarly, and also for N = 3 no 



significant difference is observed ||53|| . Another tree-level improved action is the perfect 
action ||. Also in this case numerical simulations show a better behaviour, see Fig. ||. 
Tree-level improvement seems really to work. 

Let us now compare the various improved actions in more detail. While the theories 
defined by ( J2.1| ) and ( J2.4|) behave similarly, the perfect action is vastly better and indeed it 
does not show any corrections to scaling (at the level of 1-2%) even for L = 5. A possible 
explanation of this behaviour has been suggested by Hasenfratz and Niedermayer |TU 



They show that the classically perfect action used in the simulations reported in Ref. 
|| has very small one-loop corrections so that it can be effectively considered one-loop 
improved. As an indication of how much one-loop improved an action is, they suggested 
considering the mass gap of states with spatial momentum p in perturbation theory at 
one loop. The general analytic calculation of this quantity at one loop is reported in App. 
O. If we define the correlation function 



G(t,p) = j J2 ( ff <M* ■ °W e^~^ , (5.1) 

and the mass gap uj{p) 



L 

X\,X2 



W (p)=-.im!2^M, (5.2) 

|t|— »oo t 



then, we find in perturbation theory 

1 1 



u{p)L = pL + ^+ jf{pL, p) + 0(/T 2 , L- 4 ). (5.3) 



For the various action we obtain: 



1 7T 

standard Q(pL, (3) = - — (pL) 3 + — pL; (5.4) 

(vL) 3 

on-shell Q(pL, /3) = [0.01604(AT - 1) + 0.01321] ; (5.5) 



(vL) 3 

Symanzik Q(pL, (3) = 0-01237; (5.6) 

P 

perf. k = 2 Q(pL, (3) = [a pf (N - 1) - 0.0004] . (5.7) 
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We have been unable to estimate the constant a p f precisely, but our results indicate 
\ a pf \ ~ 0.0005. From the estimates of Ref. [[HJ we would obtain the bound \a p f \ < 10~ 4 . 



As expected, for tree-level improved actions, there are no (3° corrections. Considering 
now the 0{f3~ 1 ) terms, we see that the perfect action has much smaller corrections com- 
pared to the Symanzik and the on-shell one. Therefore, if the good scaling behaviour of 
the perfect action is related to the fact that it is effectively one-loop improved — in other 
words, if higher order corrections play little role — our results could have been expected 
on the basis of the results for the mass gap we reported above. A check of this argument 
consists in verifying that an action that is exactly one-loop improved and is sufficiently lo- 
cal has the same good scaling behaviour of the perfect action. A simulation is in progress. 
Of course, at the end of this discussion, the most important question concerning the effec- 
tiveness of perfect actions remains unanswered: had we to expect a priori that the perfect 
action has very small one-loop corrections? For k = +00 the action is exactly one-loop 
improved, but the couplings are long-range so that this fact is of no practical interest. 
The interesting question is what happens for those values of k that correspond to suffi- 
ciently local actions. To answer it we have computed the four-spin coupling of the perfect 
action for various values of k and correspondingly the correction Q(pL,/3). The results 
for the perturbative coefficients (see Eq. ( [B . 3 1| ) ) are shown in Fig. where, for each 
value of k, we report two points corresponding to two different truncations of the action, 
corresponding to including respectively 51 and 771 different couplings in the four-spin 
term of the action (see App. |Cj for details). One immediately sees that the corrections 
increase strongly as k — > 0, and, for instance for k = 0.75, the corrections are of the same 
order of those of the action ( |2.1| ). In other words, the renormalizat ion-group procedure 
which is the basis of the perfect action approach does not provide naturally actions that 
are "numerically" improved: there are classically perfect actions that are relatively local 
and that behave no better than (|2.1|) . In conclusion we do not see any theoretical reason 
for the exceptionally good behaviour found in Ref. |§ . Indeed it is not a a priori obvious 
- in the sense that it is not built in the approach — that there exists a value of k such 
that the action is local and the corrections are small. 
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A Finite-size scaling in the large- N limit 
A.l Definitions 

In this section we will report some definitions and expressions that will be used in our 
computation of the corrections to FSS. 
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FSS functions are expressed in terms of the so-called remnant functions [54]. We define 

1 'k - 1/2 N 

m=0 

1 



GJa) 



E 



n=l 
oo 



(n 2 + a 2 )fe -l/2_^ 



m 



a 2m n 2k-2m-l 



H k a) 



E 



J {n 2 + a 2 ) fc+1 /2 



(A.l) 
(A.2) 



where k > (resp. fc > 1) in the first (resp. second) case. Asymptotic expressions and 
various properties of the remnant functions are reported in Ref. [3 and in Appendix A.l 



ofRef. 48 



Our results are expressed in terms of the functions F (z; p) and Fi(z;p). They are 
defined by the asymptotic expansion of the lattice sum 



L-l T-l 



IT pi + z 2/L 2 

n x =0n y =0 



(A.3) 



for L, T — > oo with z and p = L/T fixed. Here p 2 = 4sm 2 (p x /2) + 4 sin 2 (p y /2), p x 
27in x /L, and p y = 2im y /T. We obtain 



2lT 



^—\ogL + ^F 1 (z;p)+0(L-HogL). (A.4) 



Explicit expressions for F (z; p) and F\(z; p) for arbitrary p are reported in App. B of 
Ref. [fig] . Here we report them only for the strip case, i.e. for p = oo. We have 



FoO; oo) 
Fi(z\ oo) 



1 If 1 . TT 2 „ ( Z 

2i + 27^-2 log y + Go fc 

l(L-G (-W - I Z * H (— 
6 V 12 1 V27T/ / 16 192tt 3 1 \2tt 



(A.5) 



16tt 



1, TT 2 2„ ( 



2vr/ 3 



We will also need the asymptotic expansion for z — > 0. We have 



with 



F {z- oo) = ±- + F 00 + 0(,2 2 ) 



F °° = 2~rr \ lE ~ l0g ^ + 2 l0g 2 



(A.6) 

(A.7) 
(A.8) 



where 7^ ~ 0.577215664902 is Euler constant. 



A.2 Computation of the FSS functions 

In this Appendix we will compute the FSS functions and their leading corrections in the 
large- N limit for the theory defined by the action ( [3.1|) . We will follow closely the strategy 



of Ref. [ESI to which we refer for the derivation of the basic results. 
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We will consider a square lattice of size L x T. The FSS limit corresponds to (3 — > oo, 
L, T — > oo with L/T = p and m^^L = z fixed. 

The first step is the derivation of the asymptotic expression for ujl,t- From the gap 
equations ( |3.5|) and ( |3.6|) we obtain 

1 V" K{p) 
1 lt ^w(p;u L! t) +m 1 LT 
-Ul,t = i • ( A - 9 ) 



2« 3 ' i V 



LT Z^w{p ]ujLT ) +m LT 

From this equation we immediately see that, for m^y — > 0, the r.h.s. goes to zero for any 
L,T, including L,T = oo. Therefore in the FSS limit u LjT — > 0. Now, at least for 
sufficiently small, w(p; ul,t) > everywhere since w(p) > 0; moreover w(p; lol,t) vanishes 
only for p = in the Brillouin zone. Therefore, we can use the results of Ref. |^8| and we 
can write 

n -2a 3 Si(cu L , r ) /log A 



where 



S lH = / ^-^M,, (A.12) 



(27r) 2 \ u;) p 2 

(27r) 2 uj) 

and -F (,z; p) is defined in Eq. (|A.5|) . 

Eq. ([A.lOj) shows that u>l,t goes to zero in the FSS limit logarithmically and that it 
admits an expansion in powers of 1/logL. Explicitly we have 



47ra 3 S 1 (0) 

UL,T ; j 

logL 

where 



2a 3 ~ 1 (0) + 2« 3 H 2 (0) - F (z; p) - A (0)) 



logL 



-0(log- 3 L), 
(A.13) 



-2 (w) = / (A. 14) 

Eq. (|A.1C|) defines ujl,t implicitly as a function of z, L, and p, and it should be solved 
exactly if we want to obtain the corrections to FSS up to terms of order log 2 L/L 4 . In 
the following we will obtain the corrections to FSS parametrized in terms of z and ul.t] 
it is understood that Ul,t, for each z, L and p, is the solution of Eq. ( |A.10| ). 



Let us now go back discussing the validity of Eq. ([A.10Q . As already mentioned above, 
it is valid only for those values of lol,t such that w(p; uj) > 0. A second limitation derives 
from the asymptotic expansion of the denominator in the r.h.s. Indeed, if w(p; uj) > 0, 
the sum ^2 p (w(p; uj) + m^ j)" 1 is always positive. On the other hand the asymptotic 
expansion is negative for large values of z. In particular, for each L and co>l,t there is a 
unique value z c (L,ujl,t) such that 

— logL + F (z c (L,u L>T ); p) + A (u l ,t) = 0. (A.15) 

2,71 
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It is easy to check, using the results of Ref. p5| , that the asymptotic expansion makes 
sense only for z <C z c {L,uj L)T ) and that z c (L, ujl,t) behaves as L for L — > oo. In other 
words the asymptotic expansion of the lattice sum in powers of 1/L 2 is not uniform in z: 
the error increases as z — > oo. This fact is not unexpected. Indeed since for L — > oo we 
find z c (L, ujl,t) ~ L, our discussion tells us that the expansion is valid only for L ^> z, i.e. 
for rai^T *C 1, that is when the correlation length is much larger than the lattice spacing. 
In Fig. [| we report 10l,t for two different values of L and p = oo. The non- uniformity is 
clearly evident: the expansion worsens as z —>■ oo at L fixed. 

A second important fact that is evident in Eq. (|A.10|) is that the perturbative (PT) 
limit z — > at fixed L followed by the limit L — > oo does not commute with the FSS 
limit. 

Let us consider for simplicity the strip case p = oo and use the asymptotic expansion 
( A.7|) . Then, in the perturbative limit, we obtain 



-4a 3 Si(0),2 



7T 



(logL + O +0(z 3 log 2 L), 



with 



c w = 4vra 3 [S^O) - S 2 (0)] + 27rF 00 + 27rA (0). 



(A.16) 



(A.17) 



In the PT limit ul,oo has an expansion in powers of z such that the coefficient of z n is of 
order (logL) n_1 . This expansion is clearly different from the expansion ( JA.13 ), that, for 
z — > 0, becomes: 



-47ra 3 Si(0) 
logL 



7T 



z log L 



[l + 0{z)) + 0{\og- 2 L) 



(A.18) 



In the FSS limit we have a different expansion in powers of 1/ logL with coefficients that 
diverge for z — > 0. 

The next step in the derivation consists in obtaining the relation between ujl,t and 
Uoo. Starting from the second gap equation (|3.6|) we have 



1 + uj l ,t)^l,t f d 2 p 



k( P ) 



l + ^oo)^ J (27r) 2 w(p; ^oo) + 7774 LT ^ w(p;u LtT ) + m\ T 



-E 

TT 



k{p) 



(A.19) 



In the FSS limit, assuming always ujl,t and small enough so that w(p;u) > 0, using 
the results of Ref. |48|] , it is easy to compute 



-E 

TT ^ 



k{p) 



LT ^ w(p; Wi t) + m 



2 



L,T 



— logL + F (z; p) + H 3 (u; L>T ; 

Z7T 



d 2 p 



#(p) 



? 2 



--log^ 

4tt S 32 



+ S 3 (cu c 



(2tt) 2 w(p; Woo) + ™ 2 
with corrections of order m^ T logmL,T and logm^ respectively. Here 



(A.20) 
(A.21) 



d 2 P ( k(p) 1 



f2vr) 



tu(p;u;) 



(A.22) 
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Using Eq. ( A. 19 ) we obtain 



^oo = ^l.t ( 1 + A^z^l^t; p; L)— + 0(log 2 L/L 4 ; 



(A.23) 



where 



l + cu 



(1 + tj)Ei(u}) + uE 2 (uj) 



{z 2 - 32e-^ F ^) (±- log L + F (z; p) + H 3 (a;)) . (A.24) 



In the calculation of this expression we have already used the FSS result for m^/m^j, 
that will be derived in the following. The function A w (z, lo\ p; L) scales as log L for L — > oo 
and it has an expansion in powers of 1/ logL. Explicitly 



A w {z,u L)T ]p;L) = ^ 



k=-l 



\oR h L 



(A.25) 



The first coefficients are 



>;p) 



2ttH 1 (0) 



(z 2 - 32e~ 47rF " {z ^) , 



(A.26) 



M*5P) = ^^-32e^ F ^)[2a 3 E 1 (0)+F (z-,p)+E 3 (0)]. (A.27) 

Finally let us compute the corrections to FSS for the ratio m 2 O0 /m 2 LT . Using Eq. ( |3.5| ) we 
obtain 



1 



d 2 p 



1 + ujt.tLT ^ 



1 



1 + ^oo J (2tt) 2 w{p\ uJoo) +m 2 00 1 + uj L:T 
Using the results of Ref. (48| and Eq. (|A.23|) , we obtain 



w(p;u; LiT ) + m 2 LT 



(A.28) 



m 



L,T 



f m (z; p) 1 + A m (z, u LjT ; p;L)— + 0(\og 2 L/L 



where 
and 

A m (z,u LiT ;p; L) 



fm{z\p) 



32 



V- 



-4nF (z;p) 



(A.29) 



(A.30) 



-(12«! + 16d 2 - 1) (32e- 47rFo(z;p) - z 2 ) logL + 



16^(126;! + 16d 2 - l)F (z;p)e- 47TFo{z ' p) 

-4(8dx + 8d 2 - l) e - 4 ^ ^ - 4tt (jFi(z, u LiT ; p) + 32 e - 47rFo ^A 1 K iT )) 



1 + ^LT V27T 



log L + F (z;p) + E 3 (oo LiT ) A (z,ul,t',P',L). (A. 31) 
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Here 

Fi(z,u;p) = (l-12a 1 )F 1 (z;p) + z 2 (^-A 1 (u)\+2a 2 z 2 F (z;p) + ^-z 3 ^(z;p), 

' (A.32) 

Ai( ^ = /|^(w-w + w^^^ (A - 33) 

The variables &i and a 2 are functions of uj defined by the asymptotic expansion for p — ► 

«>(p; w) = p 2 + a x Y,Pt + &2{p 2 f + 0(p 6 ). (A.34) 

The function A m (z, uj; p; L) scales as logL for L — > oo and it has an expansion in powers 
of 1/ log L. Explicitly 

A m (z,uj L , T ;p;L) = E^T^- (A.35) 

The first coefficients are 

£m,-i(*; P) = ™(12<*i + 16a 2 - 16« 3 - 1) (z 2 - 32e- 4 ^ Fo(2;p) ) (A.36) 

9 

5 m ,o(z;p) = -— (l-12a x ) + — (l-12ai-12a 2 + 16a 3 ) + 0(z 2 ). (A.37) 

io 4 

One immediately sees that if the action is tree-level improved 8 m ^\(z; p) = 0. The next 
coefficient S mi0 (z; p) instead does not vanish in general. In Fig. |^ we report the function 
A m (z, ujl,t; P\ L) for p = oo together with the exact deviations from finite-size scaling. 
Notice that for L — 10 the higher-order terms in 1/L behaving as logL/L 4 still play an 
important role for z ~ 2 — 10, while they are negligible for L = 30. 

Since our class of actions for 0:3 7^ is only on-shell improved, we expect to find an 
improved behaviour only in on-shell quantities. On the strip (p = 00) let us determine 
the FSS behaviour of the mass gap p(L). In the FSS limit with Lp(L) = x fixed we have 

//! Xi ~ -/„(*) (l + A,(x,u L ^,L)j- + 0(L-Hog 2 L)] (A.38) 



where / M (x) = f m {x) and 

A^x.uj; L) = A m (x,w;oo;L) + ^(12« 1 + 12a 2 -l)e- 4 " Fo(a:;oo) + 

+(12ax + 12a 2 - 1)— (A.39) 

For L — > 00, Ap(x, uj; L) has an expansion in powers of 1/logL: 

A»(x,uj L , T ;L)=f^ 6 -^. (A.40) 
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It is easy to see that <5 Mi _i(x) = 5 m _i(x;oo). Therefore the leading term vanishes for 
tree-level improved actions. If a\ = ■h and a 2 = «3 we have 



SUx) = { 8na 2 [A (0) - 2a 2 S^O) - 2 S 2 (0) - 2 ! 



- 47T 



Ai(0) 



96tt 



1 



3(A--)+4(/3 2 -a 2 ) 



— 47ra 2 
= -47ra 3 :Ei(0) 
where fti and /3 2 are defined from 



Sx(0) J (32e- 47rFo ^ p) - x 2 ) , (A.41) 



y - ^) - x U - 1 - 327r 2 a 3 ) + 0(x 2 ) 



tf(p)=^ + Ap 4 + &(pT + 0(p b ). 



(A.42) 



(A.43) 



From Eq. ( |A.41| ) we immediately see that it is possible to choose w(p) and K{p) in such 
a way that 8^q{x) = 0. Actions satisfying this condition are one-loop on-shell improved. 
On the other hand, from Eq. (|A.42|) , one immediately convinces oneself that 8^ t i(x) never 
vanishes unless «3 = 0. In other words, actions of the form ( |3.1| ) with ^ cannot 
be two-loop improved. This is not unexpected since six-spin couplings are needed for 
two- loop improvement. 



B Perturbative computation of the mass gap for spa- 
tial momentum p ^ 

In this Appendix we will compute the mass gap for states with spatial momentum p for 
the most general action with two-spin and four-spin couplings: 

S(cr) = - w xy (T a, • <r y + c xlX2 . X3Xi {a xi ■ a X2 — \){a X3 ■ cr X4 — 1) . (B.l) 



x,y 



X1X2X3X4 



Notice that the couplings c XlX2 - X3X4 are well defined only for X\ 7^ x 2 and X3 7^ X4. This 
action reduces to Q3.1] ) if 



(B.2) 



Cx\X2',X3X4 -K XlX2 K X3X4 (8 XlX3 ~\~ 8 X2X3 ~\- S X ±X4, 8 X2 

o 

We introduce the Fourier transform 



c(p,q,k) — ^ c XlX2 - tX3X4 e 2 ^ 



+X2— X3— X4) iq(xi-X2) ik{xz— X4) 



(B.3) 



X2X3X4 



and define 



C(p;q,k) = c(p',q,k)-c(p;p/2,k)-c(p;q,p/2)+c(p;p/2,p/2), (B.4) 

p — q p — q y 



Sip, q) = c [p + q 



I, r, 5 



(B.5) 
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For p, q, k — > 0, keeping into account the lattice symmetries, we have 



-\v [(p ■ # + (p ■ m - \i £(« + + 

+aq 2 k 2 + p(q • k) 2 + 7 £ <&% + 0(p 6 , p 4 q 2 , . . .) (B.6) 



f(p,q) = a(p-q) 2 + ^[(p-q) 2 +p 2 q 2 ] +1 J2plti + (P 6 >P 4( l 2 >--^ ( R7 ) 

where a, (3 and 7 are free parameters. If c XlX2 - X3X4 is given by Eq. flB.2|) , then (3 = 7 = 
and a = —a^/2. The general conditions that make the action (B.l) tree- level on-shell 
improved have been discussed in Ref. [10|. The two-point function is improved if a± = 
1/12, while improvement of the four-point function gives 

a = ~a 2 , /3= 7 =0. (B.8) 
Let us now compute the mass gap. Explicitly we consider 

G ^P) = lJ2 ■ trt^e**'-* (B.9) 

X\,X2 

on a strip 00 x L with periodic boundary conditions in the spatial direction. For large \t\ 
we have 

G(t;p) ~ e" 1 ' 1 ^. (B.10) 

The calculation at one loop for p = is reported in Ref. f|[55||. Here we will repeat the 
calculation for p ^ 0. At the order we are interested in we expand 

oo(p) = uj (p) + 7^1 (P) + 0((3- 2 ). (B.ll) 
The tree-level term lvq(p) is easily computed. Indeed at tree level we obtain for p ^ 

G(t;p) = — — / r — . B.12 

w(q, p) is a sum of trigonometric functions and therefore it is analytic in the whole complex 
g-plane. It is then easy to compute the integral using Cauchy's theorem. For \t\ -^oowe 
have 

/V — 1 P -^o(p)|*| 

G((;p) = — -TW' (B ' 13) 

where iuo(p) is the zero of w(q,p) with the smallest positive imaginary partQ and 

D( j, ) = - i (^) . (B.14) 



V 9q 



q=iu)Q (p) 



2 We assume here that uJoijp) is real. This is not generically true, although it is always verified for 
0. 
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For p — ► we have 



D( 



T>+ [an- -^)p 3 + 0{p 5 ' 



(B.15) 
(B.16) 



For tree-level on-shell improved actions we have a% = 1/12. Therefore for this class of 
actions 0(a 2 ) corrections vanish as expected. 

Let us now consider the one-loop correction. In this case one should pay special 
attention to the boundary conditions in the temporal direction. As usual we consider 
a lattice of size L x T with free boundary conditions in the temporal direction and we 
will use the explicit expression for the two-point function reported in |55|]. G(t;p) is then 
obtained taking the limit T — > oo. For p ^ we obtain the final expression 



uj x {pL)L 



L 



1 

2 + D 



+4(JV-1) 
L 



dq w(p + q) - w(p) 



2tt 



+1 



D( 



D(p) 

tJ2 



w(q) 
dq 1 



- 1 



C(0;p,q) 



L 



dqp 1 
27r 



27r 

(/(p,?)-/(p,o)) 



+ 



+ 



(B.17) 



where in the r.h.s. p = (ico> G°) , p) • We can then expand Eq. QB.171) in powers of 1/L for 
p — > with p~L fixed. We obtain 

1 Q(pL) 



io x {p)L 



L 2 



0(L~ 



(B.18) 



where 

n(pL) 



7T 27T 

—pL (1 - 12«! - 8a 2 - 16a - I67 - 8f3) (N - l)pL{/3 + 7) 

12 3 



- Oil 



1 

12 



(N - l){pLfE + (N- l){pLfD 



A-3B-F-2[a 1 - 
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°-4 C 



The constants A, B, C, D, E, F, and G are explicitly given by 



A = 

B = 
C = 



1 

48 
1 

24 

1 

2 



d 2 q 



(2n) 2 w(q) 



-T 



d 2 q 



(2n) 2 w(q) 



d 2 q 



(2n) 2 w(q) 



-T 



dq 2 dq 2 

d 2 w . , 
— (?) 



dq 2 dq 2 
9q 2 u 



(B.19) 

(B.20) 

(B.21) 
(B.22) 
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dp 2 
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p=0 



d 2 C 
dpj 



(0;p,q) 



p=0 



(2ir) 2 w(q) 
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For the action Q3.1|) , these expressions simplify becoming 
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(2vr) 2 



d 2 K 



d 2 K 



2 (?) - 





4£ 



(27r) 2 w(g) 1 4 [ dql 
dKd 3 K OK 



dq\ 
dK d 3 K 



1 2 


" d 2 K ' 




_dq dqi_ 



G 



as 
' 4 



d 2 q 1 
(2vr) 2 w(g) 



dK 
dq 



(?) + 



9qi 



(?) 



with /3i defined in Eq. ( [031) . 

If the action is tree-level on-shell improved Eq. (|B.19|) reduces to 

Q(pL) = (pLf [(N — 1)D + A — SB — F] . 



(B.23) 
(B.24) 
(B.25) 
(B.26) 



(B.27) 
(B.28) 



(B.29) 
(B.30) 



fB.31) 



Notice that if we require the action ( |B.19|) to be tree-level improved to order 0(a A ) we 
obtain the additional condition f3\ = 1/12. In this case we also have D = 0. 

Computing numerically the various integrals, we obtain for the various actions we have 
introduced in the text: 



1. Standard action (2.3) 



7T 



2. Symanzik action 

3. Diagonal action (|2.5|) 



n(pL) = -0.01237(pL) 3 - 



7T 



Q(pL) = ^-pL - 0.00315(pL) 3 ; 
18 



(B.32) 
(B.33) 
(B.34) 
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4. On-shell action O 



n(pL) = -0.01321(pL) 3 - 0.01604(iV - l)(pLy. (B.35) 

In particular for N = 3 and pL = In we obtain in the four cases: f2(27r) = 1.64493, 
—3.0698, 0.3163, —11.235. For the first two cases we can compare with the numerical 
results of Ref. |TD[ finding good agreement. 

We have also computed Q(pL) for the one-loop Symanzik action of Ref. 0. In this case 
we should add to Eq. ( |B.19|) the tree-level contribution of the 1/(3 corrections appearing 
in the action. As expected the sum of the two terms vanishes. We should note that 
this cancellation does not happen for the Symanzik action that has been used in the 
simulations [[5(3] and that is therefore clearly incorrect. 



C Mass gap for the perfect action 

In this Appendix we want to compute the mass gap for p ^ at one-loop for a generic 
perfect action, generalizing the results of Ref. [TO]. The first step consists in deriving the 
perfect action for k ^ 2. For our calculation we will be only interested in the two-spin 
and in the four-spin couplings. We will therefore consider the generic model ( |B.1|) . With 
the standard parametrization er x = (<p x , yl — </>^.) we obtain 

{<t> Xl ■ <t>x 2 ){<t>x 3 ■ <t>x 4 ) + ■ ■ ■ (C.l) 

x,y X1X2X3X4 

where 

Ca;iX2;a;3X4 Cx±x 2 ;X3X4 &xix 2 ^ ^-x\z\x 3 X4 &x 3 X4 ^ ^-x\X2\x 3 z 

z z 

~^~^x\X2 ^x 3 X4 ^ ^ C-x\z\X3z' ~^*U XlX3 5 XlX2 5 X3XA . (C'2) 

zz' 

It is easy to see that c XlX2 - X3Xi satisfies the constraint 

^ ^ C-xiz\X2X3 g tU XlX2 5 X2X3 , (C.3) 
z 

that is related to the 0(iV)-invariance of the theory. In other words, an action of the form 
( |C . 1| ) is 0(iV)-invariant up to terms of order </> 6 if and only if Eq. ( |C.3| ) is satisfied. 

The couplings w xy and c xlX2 . )X3Xi are obtained requiring the action to be a fixed point 
of a family of renormalization-group (RG) transformations || labelled by a parameter k. 
One finds that w xy is the perfect laplacian defined in Eq. ( |2.7| ). The four-spin coupling is 
the fixed point of the equation 



— / 

C z 1 z 2 ;z 3 Z4 



^ ^ c xiX2\x 3 X4T XlZl T X2Z2 T X3Zs T XaZa + b Zl Z2;z s z 4 ■> (C-4) 



X\X 2 X 3 X4 
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where the matrix T xz is defined in Eq. (28) of Ref. || and called there M{n, ne). It 
satisfies the properties 

]TT XZ = 4, (C.5) 

X 



E v, 



(C.7) 

0, (C.8) 



x 2 ;x 3 X4 



X2X3X4 



where n is an arbitrary lattice vector. It is also important to notice that T xz is strongly 
peaked around (x — 2z) fl e {0, 1} for the range of k considered in numerical computations. 

The first step in trying to solve Eq. ( |C.4| ) (more precisely the equation c' = c) is 
to decide whether to work in real space or in Fourier space. Since Eq. (|C.4[) , once 
written in Fourier space, does not map continuous functions into continuous functions, 
the more convenient choice is the first one. Moreover by working in real space one can 
take advantage of locality. 

In order to solve Eq. ( |C.4|) one must face various technical problems. First of all 
the equation involves an infinite number of couplings. However, as noticed by Hasenfratz 
and Niedermayer ||, the relevant couplings are of short-range type unless k is very small 
or very large. One can thus hope to obtain reasonable approximations if one considers 
truncations that involve couplings among the spins of nearby points. We have therefore 
considered only couplings that are defined in an I x I square (a more precise definition 
will be given below). Notice that, even for small /, the number of couplings turns out 
to be quite large. Indeed it increases roughly as (2/ + l) 6 . In particular the number of 
inequivalent couplings is 51, 771, and 5329 respectively for I = 1,2,3. 

To define the domain precisely, let us begin by noticing that the couplings c xlX2 . X3X4 
should be invariant under the group of lattice symmetries (including of course the trans- 
lations) which acts on the four points [x] = (21, x 2 ; x 3 , x 4 ) and under the group of permu- 
tations generated by 

(x 1 ,x 2 ;x 3 ,x 4 ) -> (x 2 ,xi;x 3 ,x 4 ) (x 1 ,x 2 ;x 3 ,x 4 ) -»• {x 3 ,x± ;x 1 ,x 2 ). (C.9) 

To each set of points [x] we associate the quantity \[x]\ = max fti {i ij(i } and a new set of 
lattice points [d ] = (c?i, d 3 , d±) that are equivalent to [x] under the above symmetries 
and that satisfy the following two properties: (i) d\ = (0,0); (ii) \[d}\ < \[y]\ for all [y] 
equivalent to [d] - - and therefore to [x] - - such that y\ = (0,0). The truncation to the 
I x / square is defined by keeping all couplings c XlX2 - x . iXi such that \[d] \ < I. 

Because of this truncation, solving Eq. ( |C.4| ) becomes a standard linear algebra prob- 
lem. Eq. (p.4|) can be rewritten as 



Ci = J2^j + ^ (CIO) 

where the indices i,j run over the set A' of inequivalent (with respect to the above 
symmetries) couplings q that are defined in an I x I square. 
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Notice a further difficulty in Eq. ( U.4 ): in principle we should sum over Xi ranging on 



all Z 2 . However thanks to the observation which follows Eq. ( |C.8|) , the contributions to 
the sum are very small unless < (2/ + 1). Practically we found < {21 + 4) to be 
enough at our level of precision. 

Now we have to solve the linear system ( |C.10| ). One could try to find the solution 
iterating Eq. ( |C.10|) a sufficient number of times (after all, it is an RG transformation). 
Unfortunately this approach does not work: because of Eq. ( |C.5| ), T has an eigenvalue 
equal to 256 so that the iteration of Eq. ( |C.10| ) diverges badly. This eigenvalue is associ- 
ated to the operator: 52 xlxaxaXi ((i> Xl • <t> X2 ){<t> X3 ■ (f> x J (a rather dummy one). 
A better approach is to solve the equation 

(1-T)c = 6, (C.ll) 

with a standard method (we used both LU decomposition and singular value decomposi- 
tion). 

The problem is now that T has two unit eigenvalues 3 which are associated to two marginal 
operators: 

Ot = J dx [</>■ D0] 2 , (C.12) 
2 = J dx{4> a U(t) b )(ct) a U4) b ). (C.13) 

In other words, Eq. ( |C.4| ) does not determine c XlX2 - X3X4 uniquely, but there is the freedom 
to add two additional conditions. However there are additional constraints on c xlX2 - X3Xi 
that follow from Eq. ( pT3 ) and the fact that w xy is given. For instance, using Eq. ( |C.3| ), 
and the fact that the two-spin coupling has the standard normalization E x w xy(% — V) 2 = 
—4, one can see that the four-spin couplings should satisfy the conditions 

2j c XlX2 - X3X4 [(x 3 - xif + (x 4 - xif] = -1, (C.14) 

■t'4 

^2 ^Bl*2i*3*4 t 2 ( X 3 - 2Ti) • (X 4 - X 2 )] = -1. (C.15) 



X 2 X 3 X4 



X 2 X 3 X4 



One can check that Eqs. ( |U.14j) and ( |C.15 ) are not automatically satisfied by Eq. ( C.4 ). 



Therefore they provide the additional equations necessary to have a unique solution to 
the linear system ( |C.10|) . The solution is obtained using the singular- value decomposition 



method. We have verified that the solution of the equation (|C.4|) and of the two constraints 
flCT4]) and ( P35|) satisfies Eq. 



With this procedure we have been able to compute the couplings c XlX2 . X3XA for various 
values of k and for / = 1 and 1 = 2. In order to use the results of appendix [FJ we have 
reexpressed the integrals A, B, D, and F in terms of real-space quantities: 



3 



A necessary condition for Eq. (C.ll) to have solutions is that Vi ■ b = 0, where v\ and V2 are the left 



eigenvectors of T with eigenvalue 1. We have explicitly verified this condition. 
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Table 3: Results for the integrals appearing in the calculation of the mass gap for the 
perfect action. Here k is the parameter appearing in the RG transformations and I refers 
to the domain in which the couplings are considered. 



X 

D=-T 
6 ^ 



2 2 

0Wx0^0 X l5 



(C.17) 



Cx\X2X-$xfi ''X3X4 



X2X3X4 



I J2(x 2 - xOj - 3(x 2 - Xl )l(x 2 - x x f x | , (C.18) 



X2X3X4 



F* ^ ^ ^ {CxiX2X-}X4 (Gr XlX:i G XlX4 G X2X3 -\- G X2X4 ) 
3X4 

\ - Xi)^ - 3(X 4 - Xi)l(Xi - Xi)\ J, 



(C.19) 



where G xy is the lattice propagator, defined by G xz w zy = 5 xy . The final results are 
reported in Table [| 
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Figure 1: fi(2L, /3)2L at N = oo with fi(L,[3)L = 1.0595 fixed for the various actions 
introduced in the text. "Sym", "onshell", "std", "diag" refer respectively to the actions 

(121), ro), ro), roi). 
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Figure 2: fi(2L,P)2L at N = 3 with fi(L,(3)L = 1.0595 fixed. Crosses refer to the 
standard action, circles to the on-shell tree-level improved action, and diamonds to the 
perfect action. 
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Figure 3: One-loop perturbative coefficients for the mass gap for spatial momentum p ^ 0. 
Results for various values of the renormalization-group parameter k and truncation index 
I. Squares and diamonds correspond to D for / = 1 and I = 2 respectively. Pluses and 
crosses to the combination A — 3B — F for / = 1 and / = 2 respectively. The vertical scale 
has been multiplied by 10 3 . 
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Figure 4: ujl,oo for two different values of L. The continuous line is the exact value 
obtained from the gap equations, while the dotted line is the solution of the equation 
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Figure 5: Scaling corrections for p = oo and two different values of L, L = 10 and L = 30. 



The full line is the exact expression L 2 (m^ /m| 00 x 1/ f m (z, oo) — 1), cf. Eq. ( |A.29| ), the 
dashed line is the asymptotic expression A m (z; p; L), cf. Eq. (|A.31|) . 
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